Self-organized escape of oscillator chains in nonlinear potentials 
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We present the noise free escape of a chain of linearly interacting units from a metastable state over 
a cubic on-site potential barrier. The underlying dynamics is conservative and purely deterministic. 
The mutual interplay between nonlinearity and harmonic interactions causes an initially uniform 
lattice state to become unstable, leading to an energy redistribution with strong localization. As 
a result a spontaneously emerging localized mode grows into a critical nucleus. By surpassing this 
transition state, the nonlinear chain manages a self-organized, deterministic barrier crossing. Most 
strikingly, these noise-free, collective nonlinear escape events proceed generally by far faster than 
transitions assisted by thermal noise when the ratio between the average energy supplied per unit 
in the chain and the potential barrier energy assumes small values. 

PACS numbers: 05.40.-a, 05.45.-a, 63.20.Ry 



I. INTRODUCTION 

The cornerstone work by Kramers (for a comprehen- 
sive review see Ref. has instigated an ever ongoing 
interest in the dynamics of escape processes of single par- 
ticles, of coupled degrees of freedom or of small chains of 
coupled objects out of metastable states. In undergoing 
an escape the objects considered manage to overcome an 
energetic barrier, separating the local potential minimum 
from a neighboring attracting domain. 

Stochastic — i.e., noise-assisted — escape is the predom- 
inant phenomenon being studied in statistical physics. 
Then, the system energy fluctuates while permanently 
interacting with a thermal bath. This causes dissipation 
and local energy fluctuations and the process is condi- 
tioned on the creation of a rare optimal fluctuation which 
in turn triggers an escape In other words, an opti- 
mal fluctuation transfers sufficient energy to the chain so 
that the latter statistically overcomes the energetic bot- 
tleneck. Characteristic time-scales of such a process are 
determined by the calculation of corresponding rates of 
escape out of the domain of attraction. In this realm, 
many extensions of Kramers escape theory and of first 
passage time problems in over- and underdamped ver- 
sions have been widely investigated P, [?] . Early gener- 
alizations to multi-dimensional systems date back to the 
late 1960's Q. This method, nowadays, is commonly 
utilized in biophysical contexts and for a great many ap- 
plications occurring in physics and chemistry 

With this work we present a different scenario of a pos- 
sible exit from a metastable domain of attraction which 
has recently been put forward in Ref. [12| . The model we 
shall study is a purely deterministic dynamics of a linearly 
coupled chain of nonlinear oscillators. Put differently, no 
additional coupling to a thermal bath promotes the es- 
cape. Henceforth, dissipation vanishes as well within this 
set up. The underlying mechanism to create an escape 
is caused solely by the strongly nonlinear Hamiltonian 
deterministic dynamics. 

We explore macroscopic discrete, coupled nonlinear os- 
cillator chains with up to 1000 links. These may appear 



as realistic models in mechanical and electrical systems, 
in various biophysical contexts, in neuroscience, or in net- 
works of coupled superconductors, to name but a few 
A deterministic escape in the absence of noise 
is particularly important in the case of low temperatures 
when activated escape becomes far too slow. Also the 
case of many coupled nonlinear units in the presence of 
non-thermal intrinsic noise that scales inversely with the 
square root of the system size then calls, in the limit of 
large system sizes, for a deterministic nonlinear escape 
scenario. 

In the nonlinear regime an initially, almost homoge- 
neous chain is able to generate spontaneous structural 
modulations. This process proceeds in a self-organized 
manner. More specifically, due to the modulational insta- 
bility, unstable growing nonlinear modes give rise to the 
formation of coherent structures [181 ] such that the origi- 
nally uniformly distributed energy becomes concentrated 
to a few degrees of freedom. With regard to nonlinear 
localization phenomena intrinsic localized modes or dis- 
crete breathers, such as time-periodic and spatially local- 
ized solutions of nonlinear lattice systems, have turned 
out to present the archetype of localized excitations in 
numerous physical situations p^|-(23j. 

An escape is related to a crossing of a saddle point 
in configuration space, corresponding to a bottleneck [l[ 
or a transition state. The latter is associated with an 
activation energy E act to be concentrated in the criti- 
cal localized mode. We will show that the critical local- 
ized mode can be reached in the microcanonical situation 
spontaneously Thus, we encounter a self-organized 
creation of the transition state which is in clear contrast 
to noise activated escape. In particular, we demonstrate 
that intrinsic nonlinear effects on a long discrete chain of 
N units induce a transition over an energetic barrier by 
enhancing one, or several localized breather states fl5j |- 
[l7| . Due to this mechanism the initially almost uni- 
formly distributed energy is dynamically concentrated by 
use of an internal redistribution; no assistance of energy 
exchange with a thermal bath is thus needed. We show as 
well that the nonlinear mechanism of energy localization 
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may promote a faster escape dynamics as compared to 
the noise-assisted situation where the system experiences 
an enduring stochastic forcing. 

The paper is organized as follows: In the next section 
we introduce the model of the coupled oscillator chain 
and discuss in Sec. lIIII the modulational instability as the 
localization mechanism. In Sec. IIVI we proceed by focus- 
ing our interest on the low-energy modes corresponding 
to nearly equilibrium states of the lattice chain. The 
properties of localization induced by the dynamical for- 
mation of breather arrays are explored. Concerning the 
escape itself, special attention is paid to the passage of 
lattice states through a critical localized mode (transition 
state) in Sec. [V] Subsequently, in Sec. IVHI we demon- 
strate that the rate of escape may be crucially affected 
by the coupling strength. In Sec. |VI] the escape rates 
obtained under microcanonical conditions are compared 
with those found for thermally activated barrier cross- 
ings. In this context we assume not only flat-state initial 
preparations of the microcanonical system but also ran- 
dom chain configurations with a fairly broad distribution 
of the coordinates and/or momenta. In Sec. lVIIll we deal 
with the influence of the chain length on the escape pro- 
cess. We conclude with a summary of our results. 



II. COUPLED OSCILLATOR CHAIN MODEL 

We study a one-dimensional lattice of coupled nonlin- 
ear oscillators. Throughout the following we shall work 
with dimensionless parameters, as obtained after appro- 
priate scaling of the corresponding physical quantities. 
The coordinate q of each individual oscillator of mass 
unity evolves in a cubic, single well on-site potential of 
the form 

tf(?) = ^ 2 -f<? 3 . (1) 

This potential possesses a metastable equilibrium at 
Imin — 0, corresponding to the rest energy E min = 
and the maximum is located at q ma x = ^>o/ a with energy 
E max = AE = c^/(6a 2 ). Thus, in order for particles to 
escape from the potential well of depth AE over the en- 
ergy barrier and subsequently into the range q > q m ax a 
sufficient amount of energy has to be supplied. 

The Hamiltonian of the one-dimensional coupled non- 
linear oscillator chain is given by 

N f 2 ~| N 

H = E y + u fa) + f E - ««] 2 • ( 2 ) 

71=1 ^ ' 71 = 1 

The coordinates q n {t) quantify the displacement of the 
oscillator in the local on-site potential U at lattice site n 
(see Fig.[TJ), and p n (t) denotes the corresponding canoni- 
cally conjugate momentum. The oscillators, also referred 
to as " units" , are coupled linearly to their neighbors with 
interaction strength k. 
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Figure 1: Potential barrier landscape with a chain positioned 
at the bottom. The parameter values are a = 1 and ujq = 2. 

The equations of motion derived from the Hamiltonian 
given in Eq. ^fy then read: 

+ uj 2 q„ - aql - K [ q n+ i + q n -i - 2q n ] = . (3) 

Throughout this work we use periodic boundary condi- 
tions according to qi(t) = qpj + i(t). Note that in Eq. © 
the nonlinearity is solely contained in the local force term. 

For a setup with interacting strength k the barrier 
height can be estimated by assuming that only one unit 
of the chain is elongated, which yields a value AE K = 
(uiq + 2k) 3 I (6a 2 ). Hence, compared to the isolated unit, 
a unit coupled to its neighbors experiences an increase of 
the barrier height. 

For sufficiently small energy per unit of all chain mem- 
bers as compared to the potential barrier, a linear regime 
with a = holds true in the considered potential |T]), 
yielding oscillatory solutions in phase space such that 
the elongations are restricted to the neighborhood of the 
potential bottom. The corresponding linearized system 

Sq n + u>o5q n ~ k [ Sq n+1 + 5q n _ 1 - 2Sq n ] = , (4) 

possesses exact plane-wave solutions (phonons) 

Sq n (t) = uq exp[i(fc n — u t)] + c.c. , (5) 

The wave number k = 27tm/N, with integer m S 
[—N/2, N/2], and the frequency u> are related by the dis- 
persion relation 

uj 2 = oj 2 + 4k sin 2 Q \ . (6) 

This expression determines the frequency of linear oscil- 
lations in the phonon band with k € [— 7r, 7r]. 

The superposition of phonon modes causes oscillatory 
states wherein distinguished units may temporarily ac- 
cumulate energies that are comparable to the barrier en- 
ergy. However, in a harmonic potential with a = these 
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states are highly unlikely. If at all, they occur at a time 
scale comparable to the Poincare recurrence time of the 
system. 

Nonetheless, utilizing nonlinear effects with a > 0, 
an initial state near the metastable minimum is struc- 
turally unstable which mobilizes structural transition of 
the chain such that a transition state is adopted. This 
mechanism will be elaborated on in the next section. 



III. MODULATIONAL INSTABILITY 

It is well established that the formation of localized ex- 
citations in nonlinear systems can be caused by a modu- 
lational instability (24|-[23]. This mechanism initiates an 
instability of an initial plane-wave when small perturba- 
tions of non-vanishing wave numbers are imposed. The 
instability, giving rise to an exponential growth of the 
perturbations, destroys the initial wave at a critical wave 
number, so that a localized hump is formed. 

To analyze the nonlinear character of the solutions of 
Eq. |3]) a nonlinear discrete Schrodinger equation for the 
slowly varying envelope solution, u n (t), has been derived 

in am 



2iu) u n + k [u n+ i + u n -i - 2u n ] + 7 \u n f 



0, (7) 



with the nonlinearity parameter 7 = 10<2 2 /3wq. The sta- 
bility of a plane- wave solution of Eq. {?]) of the form 



u n (t) — uq exp(i6 n ) + c.c. 



(8) 



with 9 n = kn — uot can be investigated in the weakly 
nonlinear regime by assuming small perturbations of the 
amplitude uq and phase 9 n that have the form of sinu- 
soidal modulations with wave number Q and frequency 
il. One then finds for the perturbational wave the fol- 
lowing dispersion relation [28[ , [2^] : 



[Q - 2k sin(Q) sin(fc)] 
2k 



cos(fc) 



2k sin 2 ( — ) cos(fc) — juq 



(9) 



Stability of the perturbations necessitates that f2 is real. 
Conversely, if the right hand side in Eq. @ is negative 
the perturbation grows exponentially with a rate 
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Figure 2: Growth rate in dependence of the wave number Q 
for different coupling strengths «, as indicated on the graphs. 



The parameter values are k 
N = 100. 



0, Mo = 0.2, a = 1, o>o = 2, and 



values of the coupling strength for a mode with uq = 0.2, 
k = and N = 100. The inequality 
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puts a constraint on the allowed wave numbers. For rel- 
atively small coupling strength k = 0.02 the whole range 
of wave numbers \Q\ ^ it is responsible for the modula- 
tional instability, albeit with different weights. Enlarging 
k not only increasingly shifts the cut-off for the allowed 
wave numbers towards Q = but in addition makes the 
instability bands also narrower. In other words the mod- 
ulational instability becomes more mode-selective. Nev- 
ertheless, the maximum of the growth rate 

1 
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is unaffected by alterations of /?, while its position 



= 2 sin" 



'r4 

4k 



(12) 



(13) 



moves closer to zero with increasing coupling strength k. 

The way the growth rates with corresponding weights 
for perturbations at different wave numbers Q contribute 
to a mean growth rate is determined by the quantity T, 
reading: 



1/2 



(10) 

Notably, this modulational instability is possible only 
in the range of carrier wave numbers k € [0, 7r/2). Thus, 
patterns of short wave length are insensitive with respect 
to modulations. 

In the following we focus our interest on the k = 
mode. In Fig. [2] we depict the growth rate for different 



9 N/2 / 

r=|£r(Q = |,, fc = o 

71 = V 



(14) 



This quantifier is depicted in Fig. O as a function of the 
coupling strength k. The maximum around k ~ 0.023 
suggests that a sizable divergence of the perturbations is 
induced. To the left of its maximum the mean growth 
rate drops drastically while on its right the decrease con- 
siderably weakens. 
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Figure 3: Mean growth rate versus the coupling strength. The 
parameter values used are the same as in Fig. [2] 



IV. ENERGY SHARING AND FORMATION OF 
ARRAYS OF BREATHERS 

In order to enhance the energy localization in the dy- 
namics of Eq. ([3]) we propose the following scenario: An 
amount of energy Eq = E tota i/N is applied per unit 
which allows the activation of nonlinear, self-organized 
excitations of the chain. 

Thus, the chain possesses a total energy E to tai = NEq. 
For an escape to take place we have: E tota i > E act > AE. 
This inequality conveys the fact that more than just one 
unit governs the escape mechanism. The initial energy 
Eq is supplied as follows: (i) First, the whole chain is 
elongated homogeneously along a fixed position qo near 
the bottom of the well. [Notice that from Eq. ((5|) it fol- 
lows that this corresponds to a flat mode being equivalent 
to a k = plane- wave solution and thus qo = 2uo.] (ii) 
Then, the position of all units and their momenta are iso- 
energetically randomized while keeping the total energy 
a constant — i.e. E to tai — NEq — const. 

The random position values are chosen from a bounded 
interval |<? n (0)— go I < Aq and, likewise, the random initial 
momenta, |p n (0) — po\ < Ap. From Eq. ^ for a plane- 
wave solution with wave number k — one deduces that 
Po = 0. The whole chain is thus initialized close to an 
almost homogeneous state, but yet sufficiently displaced 
(Aq ^ 0) in order to generate nonvanishing interactions, 
enabling the exchange of energy among the coupled units. 

As the role of the deviation of the initial conditions 
from a completely homogeneous state for the instigation 
of the energy exchange process is concerned we observe 
that the attainment of an array of humps proceeds the 
faster the larger is the width Ap and/or Aq. More pre- 
cisely, due to the emergence of a modulational instabil- 
ity a pattern evolves in the course of time (of the order 
of t ~ 5 x 10 2 ) where for some lattice sites the ampli- 
tudes grow considerably whereas they remain relatively 
small in the adjacent regions. This feature is illustrated 



Figure 4: Spatio-temporal evolution of the energy distri- 
bution E n (t). Initially, the coordinates are uniformly dis- 
tributed (k — 0) within the interval \q n (0) — qo\ < Ag with 
an average at qo = 0.3 and width Aq — 0.01 and zero 
momenta — i.e., po = Ap = 0. This yields a total energy 
Etotai — 8.1 = 6.075A_E. The parameter values are given by 
a = 1, = 2, N = 100, and k = 0.3. 



in Fig. [4l depicting the spatio-temporal evolution of the 
site-energy: 



E n = § + U(q n ) + 1 



( q n+ i ~ q n f + ( q n -x - q n f 



(15) 

The breather states possessing a relatively high energy 
occur spontaneously at an average distance of the in- 
verse wave numbers Q^ ax , corresponding to the max- 



imal growth rate T max in (|12|) . Upon moving, these 
breathers tend to collide inelastically with others. In 
fact, various breathers merge to form larger-amplitude 
breathers, proceeding preferably such that the larger- 
amplitude breathers grow at the expense of the smaller 
ones. As a result, a certain amount of the total energy 
becomes strongly concentrated within confined regions of 
the chain. This localization scenario is characteristic for 
nonlinear lattice systems [2{|, [30l]-[37l|. 

For our simulations the set of coupled equations ([3]) 
has been numerically integrated by use of a fourth-order 
Runge-Kutta scheme. The accuracy of the calculation 
was checked by monitoring the conservation of the total 
energy with precision of at least 10~ 5 . The investigated 
chain consists of 100 coupled oscillators. 

To relate the energy localization with an escape over 
the barrier we note that in the beginning the aver- 
age amount of energy contained in a single unit, Eq = 
Etotai/ N, lies significantly below the barrier energy as 
expressed by the low ratio Eq/ AE = 0.06. Thus, a sin- 
gle unit must acquire the energy content of at least 16 
nearby units before it is able to overcome the barrier. 

For further illustration we depict in Fig.[5]snapshots of 
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t = 473 




Figure 5: Snapshots of the amplitudes q n (i) of a segment of 
the lattice covering ten oscillators, illustrating the formation 
of a localized structure with a subsequent barrier crossing 
event of the chain. Initial conditions are qo = 0.45, Aq = 0.01, 
and po = Ap = 0. Parameter values are cJq = 2, a = 1, 
N = 100, and k — 0.2. Top panel: snapshot taken at time 
t = 473: Emergence of localized structures. Bottom panel: 
snapshot taken at time t — 673: Depicted is a segment of 
the lattice with one unit located beyond the barrier. The 
dashed lines in the (q n — n)-plane designate the position of 
the potential's minimum and maximum. 



q n {t) at two different instants of time. In the beginning 
the energy is virtually equally shared among all units 
(not shown). After a certain time has evolved, the local 
energy accumulation is enhanced in such a manner that 
at least one of the involved units possess enough energy 
to overcome the barrier. The question then is, does such 
an escaped unit continue its excursion beyond the bar- 
rier or can it even be pulled back into the bound chain 
formation (q n < q max ) by the restoring binding forces 
exercised by its neighbors? On the other hand, the unit 
that has already escaped from the potential well might 
drag neighboring ones closer to or, in the extreme, even 
over the barrier. Thus, a concerted escape of the chain 
from the potential valley becomes plausible . 



V. TRANSITION STATE 

Whether a unit of growing amplitude can in fact escape 
from the potential well or is held back by the restoring 
forces of their neighbors depends on the corresponding 
amplitude ratio as well as on the coupling strength. The 
critical chain configuration — that is, the transition state 
separating bounded (q n < q ma x) from unbounded (q n > 
q ma x) lattice solutions — is determined by q n (t) = 0. The 
system of equations ([3]) then reduce to the stationary 
system of equations: 

dU 

- -^z- + K[q n +i + q n -i - 2g„] = 0. (16) 

oqn 

Interpreting n as a 'discrete' time, with 1 < n < N, 
equation (|16p describes the motion of a point particle in 
the inverted potential —U(q). This difference system can 
be cast in the form of a two-dimensional map by defining 
x n = q n and y n = q n -i 38], which gives 

x n +i = (u%x n - ax 2 n )/ k + 2x n - y n 

Vn+l = x n . (17) 

The fixed points of this map are found as 

xo = yo = o, x\ = yi = — ■ (18) 

a 

A linear stability analysis reveals that (xq, yo) represents 
an unstable hyperbolic equilibrium while at (xi,yi) a 
stable center is located. The map is non-integrable. The 
stable and unstable manifolds of the hyperbolic point in- 
tersect each other, yielding homoclinic crossings as il- 
lustrated in Fig. [6l The corresponding homoclinic orbit 
of the map, consisting of the points of intersections be- 
tween W s and W u , is on the lattice chain equivalent to a 
stationary localized hump solution {q^}, centered at site 
n — n c , which resembles the form of a (pointed) hairpin 
(for details concerning the relation between homoclinic 
orbits and localized lattice solutions see [38| and [3!|). In 
Fig. [7] profiles of this hairpin-like critical localized mode 
(c.l.m.), or critical nucleus, with displacements {q%} are 
depicted for several coupling strengths. We observe that 
the stronger the coupling is, the larger the maximal am- 
plitude of the humps is, q^ = q max , and the wider the 
spatial extension of the latter is. We underline that on a 
sufficiently extended lattice this c.l.m. represents a nar- 
row chain formation with its width being much smaller 
than the total chain length. Equation (TIT))) can be derived 
from an energy functional with vanishing kinetic energy, 
reading 

Eim] = E + \ fa - 9«-i] 2 ) ( 19 ) 

n 

with dE / dq n — 0. Apparently, with increasing coupling 
strength k sl larger activation energy 

Eact = E[{q*}] (20) 



n — n,; 



Figure 6: First intersections of the stable W s (thin gray line) 
and the unstable W u manifold (thick dark line) of the hy- 
perbolic point of the map in Eq. (11T|) at the map origin 
(0,0). The parameter values are ujq = 2, a = 1, and k = 1. 
The homoclinic point labeled (Qn c j <3n c +i) an< i its map iter- 
ate (<j^ c +i; Qn c +2) correspond on the consecutive lattice sites 
n = n c , Tic + 1, n c + 2 to a decaying pattern whose ampli- 
tudes with q„ > q~n + i are represented by the arrows of varying 
length in the schematic representation in the upper left corner. 
Map iterations of the homoclinic point (<Jn c +ii <7n c +2) result 
in further homoclinic points approaching asymptotically the 
hyperbolic point at the map origin. 



is required to bring the chain into its critical localized 
mode configuration. 

The activation energies for different coupling strengths 
are depicted in Fig. [8] For small coupling the maximal 
amplitude q^ nax is beyond but still close to q m ax = <^o/ a 
[the position of the maximum of the potential U (q)] while 
the remaining units practically reside at the minimum of 
the potential U(q). For larger couplings the maximal am- 
plitude q^ ax lies at larger distances beyond q ma x- Most 
importantly, for escape we have that dU(q)/dq\g = gh < 
0. 

It can be shown that the critical localized mode, being 
associated with an unstable saddle point in configuration 
space, is indeed dynamically unstable. Setting q n (t) — 
q'n + w n (t) with \w n \ <C 1 the linearized equations of 
motion are derived as 



w n {t) 



d 2 U(g n ) 

dql 

n\w n+1 (t) 



\q n =qhW n (t) 

+ U7 n _i(t) - 2w n (t)] 



(21) 



With the ansatz w n (t) = 4> n exp(^/\t) for the solution of 
(|21|) one arrives at an eigenvalue problem 

\4> n = -V n (j) n + k[4> 71+ i + 4> n -x - 2(f> n ] , (22) 



with 



V„ 



d 2 U(q n ) , 



dql 



2 l"?n 



=qh = w - 2aq 



(23) 



Figure 7: Amplitude profile of the critical chain configuration, 
being symmetrically centered at site n c , for different coupling 
strengths ft = 0.1 (circles), k = 0.3 (diamonds), and k = 1.0 
(triangles). For better illustration only we depict the cen- 
tral part of the lattice chain with sizable elongations. The 
parameter values are luq — 2 and a = 1. 




Figure 8: Activation energy E ac t as a function of the coupling 
strength k. Solid line: fitted data with slope s = 3.54 ± 0.01. 
The remaining parameter values are lo% — 2 and a = 1. 



The second-order difference equation (|22| is of the dis- 
crete stationary Schrodinger type, with a non-periodic 
potential —V ni breaking the translational invariance so 
that localized solutions exist (so called stop-gap states). 
The evolution of the two-component vector (<f) n +i, <j> n )' r 
is determined by the following Poincare map: 



M 



1 



-1 




(24) 



with on-site energy E n = (A + V n )/n + 2. The node- 
less even-parity ground state of Eq. (|22|) . with its energy 
under the lower edge of the energy band of the passing 
band states, corresponds to an orbit of the linear map M. 
being homoclinic to the hyperbolic equilibrium point at 
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Figure 9: The escape time T e " c of the individual chain units 
versus position for collective escape and periodic boundary 
conditions for a chain consisting of 100 units for one realiza- 
tion of initial conditions. The parameter values are a = 1, 
luq = 2, N — 100, and k, = 0.3. The initial conditions are 
given by qo — 0.45, Aq = 0.1, and po = Ap — 0. 

the origin (0, 0) of the map plane. For the presence of a 
hyperbolic equilibrium the following inequality has to be 
satisfied: 

Trace(TW) = E n = X + Vn + 2 > 2 , (25) 
re 

implying that A must fulfill the constraint 

A > max(-V„) = 2a max ct - u)l > 0. (26) 

n n 

With the maximal amplitude of the c.l.m. lying beyond 
the barrier — viz., max„ > ui^/a — one finds 

A > u% > . (27) 

Therefore, the ground state belongs to a positive eigen- 
value from which we deduce that perturbations of the 
corresponding solution in the time domain grow exponen- 
tially. Hence, if the kinetic energy overcomes the critical 
nucleus, the subsequent escape of its neighbors is ini- 
tiated, which progress on the chain to the left and to 
the right of the hair pin as a pro pag atin g ki nk and anti- 
kink, respectively (see Refs. 0,113 SHI)- In phase 
space the units move parallel to the unstable manifold of 
the hyperbolic equilibrium [which is related to the saddle 
point at the maximum of the potential U (q)] , realizing 
in this way an efficient lowering of the total potential en- 
ergy. Because the kinetic energy of this outward motion 
is consequently increasing, a return backwards over the 
barrier into the original well is hereby prevented. Fig. [5] 
illustrates the kink-antikink motion showing the escape 

(n) 

time Tesc of the units versus the position on the lattice. 
The escape time of a unit is defined as the moment at 
which it passes through the value q esc = 5 q max beyond 



the barrier. We remark that q esc is chosen such that 
U(q esc ) is sufficiently lowered so that the return of an 
escaped unit over the barrier into the potential well is 
practically excluded. Consecutively, all oscillators man- 
age to eventually climb over the barrier one after another 
in a relatively short time interval. (The position of the 
first escape event varies in general for random samples of 
initial conditions.) 



VI. COMPARISON WITH THERMALLY 
ACTIVATED ESCAPE 

We next compare the microcanonical escape process 
with the corresponding thermally assisted escape process 
at a temperature T (THi H 0, IS S3, S3, SI • The associ- 
ated Langevin system reads 

d q n dll dq n 
-aW + dq- n - K [qn+1 + ^ ^H^+^t) = • 

(28) 

with a friction parameter v and where £ ra (i) denotes a 
Gaussian distributed thermal, white noise of vanishing 
mean < £ n (i) >= 0, obeying the well-known fluctuation- 
dissipation relation < £n(t)£n'(t') >= 2vkBT5 n ,n'5{t — t') 

with ks denoting the Boltzmann constant. We define the 
escape time of a chain as the mean value of the escape 
times of its units (see again above) . 



Escape times 

Our results are summarized in Fig. [TU] depicting the 
mean escape time as a function of Eo/AE. The aver- 
ages were performed over 500 realizations of random ini- 
tial conditions or the noise in the microcanonical and 
Langevin systems, respectively. For the deterministic 
and conservative system Q the excitation energy Eg is 
given by the (average) initial energy content of one unit 
Eq = Efotai/N. For the simulations we varied go while 
keeping Aq — 0.01 and po — Ap — and obtained thus 
different values of E to tai ■ In the case of the Langevin sys- 
tem (|2"5)) for sufficiently low T the energy Eq is taken as 
the thermal energy fc^T. This ratio thus corresponds to 
the inverse Arrhenius factor indeed, at sufficient low 
temperature (large ratio AE/Eq) the logarithmic escape 
time follows an almost linear behavior versus the Arrhe- 
nius factor, as expected for a noise-driven escape at weak 
noise strength 

The Langevin equations were numerically integrated 
using a two-order Heun stochastic solver scheme [44| . In 
both cases there occurs a rather distinct decay of T esc 
with growing ratio Eq/AE in the low energy region. This 
effect weakens gradually upon further increasing Eg . Re- 
markably, for low Eq (indicated as the breather region 
in the plot) the escape proceeds distinctly faster for the 
noise-free case as compared with a situation of a chain 
that is coupled to a heat bath at temperature T . This 
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Figure 10: Mean escape time T e3C vs. the ratio Eo/AE 
of energy per unit Eq and the barrier energy AE for the 
noise-free case with initial k = mode (solid line) with 
Aq = 0.01 and po = Ap = (for go, see text), the case 
with strongly randomized initial conditions (diamonds) where 
-0.251 < q n (0) < 0.534, -0.251 < p„(0) < 0.534, and the 
noise-assisted escape (dashed line), respectively. The param- 
eter values are uIq = 2, o = 1, K = 0.3, and the friction 
strength is v — 0.1. The inset depicts the same data, but now 
plotted versus the effective Arrhenius factor AE/Eq. 



implies a large enhancement of the rate of escape as com- 
pared to the thermal rate. Near Eo/AE > 0.36, there 
occurs a crossover, with the mean escape time of the de- 
terministic system at even higher ratios closely following 
that of the thermal Langevin dynamics. At these values 
the escape times become comparable with the relaxation 
time which is determined by the inverse of the friction 
strength. Apparently in the region of lower ksT nonlin- 
ear excitations are damped out in the Langevin system 
at longer time scales. Hence they will not accelerate the 
escapes in the case with fluctuations and damping. 

To sharpen our finding that the escape proceeds typi- 
cally faster in the noiseless situation as compared to the 
case with a coupling to a heat bath, we investigated also 
the escape process of non-flat chain patterns starting out 
from strongly randomized initial conditions. For these ini- 
tial conditions the coordinates and momenta are chosen 
at random from fairly broad ranges —0.251 < (?„(0) < 
0.534 and -0.251 < p n (0) < 0.534. For various energy 
values the averages were performed over 100 realizations 
of initial conditions belonging to iso-energetic configura- 
tions with ratio Eq / AE each. 

The findings for the mean escape time as a function of 
the mean initial energy content of the units relative to the 
barrier height Eq/AE are included in Fig. [TU] with the 
diamond symbols. Most importantly, even for random 
initial conditions the mean escape time assumes smaller 
values in the microcanonical situation as compared to the 
Langevin dynamics. This underpins our general state- 
ment that noiseless escape indeed proceeds faster than 
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Figure 11: Resonant behavior in the mean escape time T e3C 
as a function of the coupling strength k with the remaining 
parameter values as in Fig.U The initial conditions are for the 
coordinates qo = 0.45 and Aq = 0.1 and for the momenta po = 
Ap = 0. The averages were performed over 500 realizations 
of random initial conditions. 



thermally activated escape. 

We note that the breathers present robust chain config- 
urations that are formed rather fast as compared to the 
escape time. In contrast, the forever impinging stochas- 
tic forces seemingly impede a fast growth of the criti- 
cal nucleus and may even cause a possible destruction of 
the critical chain formation, leading to re-crossings of the 
transition region, which only hampers a speedy escape. 
This inhibition for escape is most effective at small ratios 
of Eq I AE, being induced either by high barrier heights, 
or low temperatures (implying a small Eq). A determinis- 
tic scenario thus presents a more favorable route towards 
accelerated escape in situations with very weak noise or 
very large barrier heights. Having performed also simu- 
lations for more general situations (i) with nonharmonic, 
nonlinear chain interactions, (ii) in higher dimensions, 
and (iii) with differing on-site potentials, we find [4lJ 
that the phenomenon of an enhanced, noise-free escape 
remains robust in regimes of a large effective Arrhenius 
factor with the latter given by the ratio of the barrier 
height AE and the initial energy per unit, Eq. 



VII. OPTIMAL COUPLING AND RESONANCE 
STRUCTURE IN THE ESCAPE PROCESS 

Furthermore, we study the impact of the coupling 
strength k on the escape process. The results concern- 
ing the mean escape time are illustrated with Fig. [TT1 
Strikingly, the mean escape time exhibits a resonance 
structure; viz., there exists a coupling strength (n res = 
0.31) for which the escape proceeds faster than for all 
other couplings strengths. Upon lowering k < K res we 
notice a substantial rise of the escape time while for 
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K > K r es the graph exhibits only a moderately grow- 
ing slope with growing coupling strength k. In this sense 
K res = 0.31 from Fig. [IT] represents indeed the optimal 
coupling strength for which escape is achieved within 
a minimal amount of time. Finally, outside the range 
K € [0.05 , 1.5] not even the escape of a single unit has 
been observed. The reason is that the time scale for a 
pronounced formation of energy concentration, being vi- 
tal for escape (due to breather coalescence and energy 
accumulation in the critical localized mode) , exceeds the 
simulation time (taken here as t = 5000). 

A physical explanation for the appearance of a 
resonance-like structure can be given in terms of the 
different degrees of instability of the underlying motion 
facilitating the destruction of the initial flat mode by 
modulational perturbations. We recollect that with the 
variation of the perturbation strength the growth rate 
changes (cf. Sec. IIII[) from a more flat to a strongly 
curved single-peaked structure. To illustrate the im- 
pact of the growth rate on the degree of localization 
of emerging patterns we present in Fig. [T^] the energy 
distribution defined in Eq. (|15p at an early instant of 
time — namely, after the formation of the spatially lo- 
calized structure due to spontaneous modulational in- 
stability has taken place. For comparison, patterns for 
three coupling strengths are shown. In all cases a number 
of isolated localized humps are formed. The number of 
humps, Nhumps, can be attributed to the wave number 
at maximal growth rate as follows: \ pa tternN humps = N 
and Xpattem = ^/Qmax- Most importantly, the num- 
ber of humps (besides their height and width) regulates 
how the total energy is shared among them. Clearly, for 
k = 0.31 [Fig. [T2] (b)] the energy is more strongly lo- 
calized (fewer humps and of higher height) than in the 
case of k = 0.09 [Fig. [T2] (a)]. In comparison for k = 1 
[Fig. [T2] (c)] the number of humps is further diminished, 
but they are of lower height than most of the humps for 
k = 0.31. In fact, for k — 0.31 the energy contained in 
the unit at site n = 35 is close to the one of the barrier. 
Thus, a localized pattern appropriate for escape is pro- 
vided already by the mechanism of modulational insta- 
bility. In particular, no further (major) energy accumula- 
tion, which would delay the escape process considerably, 
is hence required. 

To gain further insight into the efficiency of energy 
localization it is illustrative to suppose that the whole 
lattice can be divided into a periodic array of (non- 
interacting) segments where each of them supports a sin- 
gle localized hump. The energy of one segment is deter- 
mined by 
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Figure 12: Spatially localized structure at t = 500 for three 
different coupling strengths. Initial conditions and parameter 
values as in Fig. [4] except for the coupling strengths (a) k = 
0.05, (b) k = 0.31, and (c) « = 1.00. 



E„ = 



E 



total 



E, 



total 



N h u 



nips 



N/Xpatt e 



(29) 



Defining e s = E S /AE as the ratio of the energy per seg- 
ment to the net barrier energy we obtain 



2irR 



total 



Q max NAE 



(30) 



Appropriate conditions for escape are provided when the 
energy contained in each segment, e s , is close to the ac- 
tivation energy, e ac t, measured in units of the barrier 
energy— i.e., e act = E act /AE. The efficiency of energy 
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Figure 13: The ratio \ defined in Eq. (|31[) as a function of 
the coupling strength k. The parameter values are loq — 2 
and o = l. The initial conditions are given by qo = 0.45 and 
Po = 0. 



Figure 14: First, mean, and last escape time T eS c as a function 
of the chain length, N. The parameter values are io% = 2, 
a = l, and k = 0.2. The initial conditions are given by 
qo = 0.50, Aq = 0.01, and p Q = Ap = 0. 



localization is then determined by the following ratio 

X=— ■ (31) 

For a given value of the coupling strength k the activation 
energy is known; cf. Fig. [5] Fixing the initial energy and 
using Eq. (|10[) we infer the value of Q m ax and finally using 
Eqs. (|30p and (|3"Tj) we obtain \- In Fig. [13] the ratio \ 
is plotted as a function of the coupling strength k. The 
plot exhibits a maximum at k — 0.3, which corroborates 
the finding of the resonance found for the escape versus 
coupling strength as depicted in Fig. [TTJ Moreover, the 
curvature of the graphs of Figs. [TT1 and [lUl are similar. 



In the case of Uq > 2k/7 the initial homogeneous state 
is always unstable, independent of the number of oscil- 
lators, N. However, for an initial state in the weakly 
nonlinear regime, which means itg < 2/c/7, the inequal- 
ity (|32p yields a condition for the minimal number of 
units on the chain that are necessary for modulational 
instability. On the other hand, once the conditions are 
provided that the chain be able to adopt the transition 
state, the addition of further lattice sites beyond a cer- 
tain number leaves the activation energy unaltered. This 
is due to the fact that the transition state is represented 
by the c.l.m. which is strongly localized in space with 
exponential decaying tails. 



VIII. INFLUENCE OF CHAIN LENGTH ON 
THE ESCAPE TIME 

We also study the influence of "size" — i.e., the num- 
ber of oscillators, N, on the escape process. A general 
constraint on the escape process arising from varying the 
chain length is formulated. We discuss the escape time 
statistics for chains with constant energy density and for 
chains with fixed total energy, but varying number of 
oscillators, respectively. Our studies apply to a homoge- 
neous initial state (k = mode). 

First of all, the initial homogeneous state must become 
unstable with respect to a modulational instability. With 
relation (TTTj) a constraint is established on the allowed 
wave numbers, giving rise to the modulational instability. 
They form a discrete set, and we can derive a lower bound 
for the number of oscillators, N m i n , needed for the onset 
of the modulational instability: 




Case with constant energy density 

Let us first consider chains with constant energy den- 
sity p = Etotai/N = const. One would at first glance ex- 
pect a faster escape with increasing number of oscillators 
and thus with increasing total energy. This is, however, 
not necessarily the case. For an explanation it is suitable 
to consider the limit of very long chains, N — > 1000. 

In Fig. [T2] the average time for which the first and last 
escape incidents of a unit take place is depicted versus 
varying chain length N. The averages were performed 
again over 500 realizations of random initial conditions. 
In addition, we as well depict the mean escape time of 
the chain. We set the initial coordinates around a mean 
value of go = 0.50 and spread at Aq = 0.01, yielding 
p = 0.156AS. Apparently the longer the chain is, the 
more humps (breathers) that are formed due to the mod- 
ulational instability. This offers the possibility that a 
larger number of interacting breathers contribute to an 
enhanced energy localization in a confined region of the 
chain which in turn boosts the formation of the critical 
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Figure 15: Mean escape time T ecs as a function of the number 
of oscillators, N, for different values of Etotai- For the choice 
of go, see text. Furthermore Aq = 0.01 and po = Ap — 0. 
The parameter values are u>o = 2, a = 1, ft = 0.2, and Etotai 
as depicted in the legend. 

localized mode. Hence, the time it takes for the first unit 
to escape shrinks with increasing chain length, while the 
last escape time increases due to the enlarged number of 
escaping units. Also, the mean escape time becomes in- 
sensitive to variations of the chain length for sufficiently 
large length N > 500, and thus it tends to saturate. 

Fixed total energy 

We next consider the situation when a fixed amount 
of total energy — i.e., E to tai = const — is provided to the 
system and the number of units on the lattice is var- 
ied. To obtain a certain value of Etotai upon altering the 
number of units, N, we adopted q$ appropriately while 
keeping Aq = 0.01 and pa = Ap = fixed. The max- 
imal energy content per unit is restricted to the range 
Eq < 0.5AE. The results for the mean escape time are 
depicted in Fig. 1151 Generally, we observe an increase 
of T esc with growing N. Interestingly enough, the slope 
passes through intermediate stages of sub-exponential, 
hyper-exponential and eventually approaches an expo- 
nential behavior. 



IX. SUMMARY 

In this paper we have explored the conservative and 
deterministic dynamics of a one-dimensional chain con- 
sisting of linearly coupled anharmonic oscillators that are 
placed into a cubic on-site potential. Attention has been 
paid to the collective barrier crossing of the whole chain. 
Initially the system is placed into a metastable state for 
which all units are trapped near the bottom of the poten- 
tial. An overcoming of the barrier of the whole chain is 
prevented at short initial times because of a too high net 
barrier height. Nevertheless, as we convincingly demon- 
strated, the spontaneous formation of localized modes 
upon evolving time serves to enrich energetically a seg- 
ment on the chain to such a degree that it adopts the 
transition state energy in assuming the form of a hair- 
pin. We have shown that the associated critical local- 
ized lattice state is dynamically unstable and eventually 
a barrier crossing proceeds as the propagation of a kink- 
antikink-like pair along the chain. Strikingly, there exists 
a resonant-like coupling strength k for which the escape 
time (rate) becomes minimal (maximal), cf. Fig. 1111 

In view of potential applications we note that this de- 
terministic collective escape process provides nonlinear 
systems with the unique possibility to self-promote their 
activation dynamics. Particularly, the ability to operate 
efficiently — i.e., exhibiting an enhanced collective coher- 
ent escape — although not optimally initialized (meaning 
that one starts out with a far too low energy density 
compared to the barrier height) underpins the beneficial 
use of this physical scenario. Remarkably, while at weak 
thermal noise the rate of thermal escape is exponentially 
suppressed, a deterministic nonlinear breather dynam- 
ics yields a robust critical nucleus configuration, which 
in turn causes an enhancement of the noise-free escape 
rate. Thus, the freezing out of noise may prove advanta- 
geous for transport in metastable landscapes, whenever 
the deterministic escape dynamics can be launched in a 
single shot via an initial energy supply. 
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